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CONTRACT  SUMMARY 


This  report  summarizes  work  performed  by  SAIC  on 
Contract  N0001 4-83-C-0303  for  the  Bottom  Interaction  Program 
managed  by  NORDA.  Four  tasks  were  addressed  in  the  course  of 
this  work: 

Task  1  -  Support  the  Bottom  Interaction  Program  - 
Thin  Sediment 

Principal  Investigator  -  C.W.  Spofford 

Task  2  -  Analyze  and  Interpret  Bearing  Stake  Coher¬ 
ence  Data 

Principal  Investigators  -  L.B.  Dozier 

J.S.  Hanna 

Task  3  -  Support  the  Bottom  Interaction  Program  - 
Thin  Layers 

Principal  Investigator  -  R.R.  Greene 

Task  4  -  Investigate  Spatial  Coherence  in  Shallow 
Water 

Principal  Investigator  -  L.B.  Dozier 

Work  on  each  of  these  tasks  is  summarized  in  the  following 
sections.  Separate  SAIC  reports  have  been  generated  to  docu¬ 
ment  the  results  of  Tasks  2#  3»  and  4.  Work  on  Task  1  is 
on-going  and  has  bee  .sntinued  in  FY85  under  a  new  contract. 
Two  references  generated  under  this  effort  are  included  here 
since  they  have  to  date  received  less  distribution. 


Task  1  -  Support  the  Bottom  Interaction 
Program  -  Thin  Sediment 

Summary 


The  objective  of  this  task  was  to  attempt  to  recon¬ 
cile  several  different  aspects  of  the  problem  of  bottom 
interaction  in  the  vast,  thin-sediment  regions  of  the  North¬ 
eastern  Pacific  abyssal  hills.  Data  from  these  areas 
suggested  anomalous  behavior  of  bottom-interacting  signals  in 
two  ways; 

1 .  Loss  versus  grazing  angle  at  the  acoustic  base¬ 
ment  -  approximately  6  dB  independent  of 
grazing  angle  and  frequency  (from  50  to  1500 
Hz);  and 

2.  Very  large  time  spreading  of  the  signal  -  again 
largely  independent  of  grazing  angle  and  fre¬ 
quency. 

In  earlier  work  for  the  Bottom  Interaction  Program  SAIC  had 
developed  a  composite-roughness  scattering  theory  which  takes 
as  input  the  bottom  roughness  spectrum  and  computes  the  frac¬ 
tion  of  energy  scattered  into  specular  and  non-specular 
directions  as  a  result  of  that  roughness.  In  this  task  the 
first  step  was  to  obtain  data  on  the  basement  roughness  and 
see  whether  the  observed  roughness  when  used  in  this  state- 
of-the-art  scattering  theory  could  account  for  the  large 
observed  time  spreads.  It  was  also  hoped  that  an  understand¬ 
ing  of  this  scattering  phenomenon  might  elucidate  the 
peculiar  total  loss  characteristics. 


The  work  on  this  task  has  produced  the  following; 

1.  A  paper  summarizing  the  status  of  the  bottom- 
loss  and  coherence  problem  with  emphasis  on  the 
thin  sediment  behavior;  "Status  of  the  Bottom 
Loss  Problem  Relative  to  Bottom  Bounce  Coher¬ 
ence"  by  C.W.  Spofford,  published  in  the  Pre- 
cedings  of  the  NAVSEA  Stochastic  Modeling  Work¬ 
shop/  NAVSEA/  C.W.  Spofford  and  J.M.  Haynes, 
ed'S.,  September  1983. 

2.  A  paper  on  the  statistics  of  basement  roughness 
based  on  our  analysis  of  Deep  Tow  data  provided 
by  Fred  Spiess  of  Scripps;  "Measurements  of 
Basement  Roughness  and  Sediment  Thickness 
Statistics  in  Pacific  Abyssal  Hills,"  C.W. 
Spofford,  E.S.  Holmes  and  F.N.  Spiess,  sub¬ 
mitted  to  Journal  of  Geophysical  Research 
(1985). 

Copies  of  these  papers  are  attached  for  reference.  The 
results  presented  in  the  second  paper  represent  a  fundamental 
extension  of  the  data  on  bottom  (or  sub-bottom)  roughness. 
In  summary,  the  acoustical  basement  in  this  area  northeast  of 
Hawaii  is  composed  of  a  rough  chert  layer  lying  approximately 
40  m  below  the  water-sediment  interface.  Our  spectral  esti¬ 
mates  from  Deep  Tow  blended  in  very  smoothly  with  longer- 
wavelength  estimates  of  bottom  roughness  reported  by  Bell. 
The  spectral  slope,  however,  steepened  considerably  (from 
k“2.5  to  k~3*^  or  more)  in  the  t\i  o  decades  in  wavelength 
that  our  new  results  spanned.  Any  attempts  to  extrapolate 
Bell's  data  beyond  its  limits  would  have  been  in  significant 


error. 


In  the  second  part  of  this  effort,  this  spectrum 
and  the  observed  slope  distributions  were  used  in  the 
composite-roughness  and  specular-scatter  models,  respec¬ 
tively,  developed  by  S.MC  to  predict  corresponding  time 
spreads.  Using  the  measured  basement  roughness,  neither  the 
low-frequency  diffractive  scattering  model,  nor  the  high- 
frequency  optical  scatterer  model  were  able  to  yield  time 
spreads  as  large  as  those  observed.  In  fact  no  reasonable 
estimates  for  roughness  which  might  be  even  plausible  given 
the  roughness  data  were  able  to  produce  long  enough  time 
spreads. 

Work  on  this  problem  is  contlijuiag  on  a  follow-on 
BIP  contract  on  two  fronts: 

1.  Horizontal  and  vertical  array  angle,  frequency 
and  time  spread  data  from  this  area  are  being 
analyzed  to  better  quantify  the  nature  of  the 
scatter  both  to  support  system  performance 
estimates,  and  to  provide  insight  into  the 
effective  scattering  kernel  of  this  bottom. 

2.  The  acoustical  properties  of  the  chert,  both 
its  roughness  scales  (perhaps  not  measurable 
from  Deep  Tow)  and  its  impedance  characteris¬ 
tics  are  being  investigated.  Certainly  over 
some  of  the  area  of  interest  chert  is  the 
effective  basement. 

In  other  parts,  basalt  is  nearer  the  surface,  how¬ 
ever  some  DSDP  data  suggest  that  there  may  be  chert  very  near 
this  basalt.  This  may  be  a  clue  as  to  why  scattering  in  the 
Pacific  abyssal  hills  is  so  much  stronger  than  in  the  Hid- 
Atlantic  Ridge  areas  as  reported  by  Ellinthorpe. 


Task  2  -  Analyze  and  Interpret  Bearing  Stake 
Coherence  Data 


A  report  entitled  "Bearing  Stake  Coherence 
Analysis"  (SAIC-85/1031 )  dated  February  1985  was  delivered 
under  this  task.  The  summary  and  conclusions  from  this 
report  are  as  follows. 

We  compared  results  of  a  mode-based  spatial  coher¬ 
ence  model  with  BEARING  STAKE  coherence  measuremencs  taken  in 
a  deep-water  environment. 

The  model  predicted  excellent  coherence  in  this 
situation;  such  excellent  coherence  is  indeed  seen  in  the 
best  of  the  data.  The  apparently  poor  coherence  in  some  of 
the  data  was  evidently  due  to  causes  other  than  bottom  inter¬ 
action. 
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TasK  3  -  Support  the  Bottom  Interaction 
Program  -  Thin  Layers 


A  report  entitled  "A  Random  Layer  Model  for  Reflec¬ 
tion  from  Ocean  Bottom  Sediments"  {SAIC-85/1036)  dated 
February  1985  was  delivered  under  this  task.  The  conclusions 
from  this  report  are  as  follows. 

A  procedure  for  estimating  reflection  losses  from 
layered  ocean  bottom  sediments,  such  as  those  observed  in 
abyssal  plains,  was  implemented.  Limited  experience  with  it 
indicates  that  it  can  reproduce  observed  reflection-loss 
levels  using  the  geoacoustic  properties  and  layering  statis¬ 
tics  of  ocean  bottom  sediments.  Furthermore,  it  reproduces 
the  observed  qualitative  dependence  of  reflection  loss  on 
angle  and  frequency,  including  the  negative  frequency  depend- 
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Task  4  -  Investigate  Spatial  Coherence  in  Shallow  Water 


A  report  entitled  "Spatial  Coherence  in  Shallow 
Water"  (SAIC-85/1 035)  dated  February  1985  was  delivered  under 
this  task.  The  summary  and  conclusions  from  this  report  are 
as  follows. 

We  compared  model  predictions  of  horizontal  signal 
coherence  loss  due  to  bottom  volume  scattering  with  data  from 
a  shallow  water  experiment.  Although  we  did  not  obtain  quan¬ 
titative  agreement  because  of  apparent  rough-surface  coher¬ 
ence  loss  included  in  the  data,  we  were  able  to  identify 
mode-stripping  over  the  same-range  intervals  in  both  the 
model  and  the  data.  In  addition,  we  developed  a  rule  of 
thumb  that,  for  this  experiment  at  least,  the  optimal  fre¬ 
quency  for  transmission  was  also  relatively  optimal  for 
coherence. 
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ABSTRACT 


This  paper  attempts  to  summarize  our  current  under¬ 
standing  of  low-frequency  bottom  loss  in  order  to  set  the 
stage  for  addressing  the  problems  of  bottom- bounce  coherence. 
The  bottom-loss  problem  in  areas  where  sediments  are  thick 
seems  to  be  well  in  hand.  This  understanding  is  demonstrated 
through  comparisons  of  measured  loss  with  simulated  loss 
using  a  relatively  simple  geo-acoustic  model.  In  thin  sedi¬ 
ment  areas  the  loss  is  modelable,  but  with  somewhat  lower 
accuracy  because  large-angle  scattering  is  so  dominant.  The 
data  collected  in  these  areas  suggest  three  basic  coherence 
problems: 

1.  Effects  associated  with  fine,  near-surface 
layering  in  the  sediments, 

2.  Coherence  losses  due  to  random  inhomogeneities 
in  the  sediment  volume  itself,  and 

3.  Large-angle  scattering  and  associated  coherence 
problems  when  the  sound  interacts  with  the 
rough,  hard  basement  in  thin-sediment  areas. 

Work  being  performed  at  SAI  in  all  three  areas  is  summarized. 
(A  subsequent  paper  by  Dozier  in  these  proceedings  discusses 
the  second  problem  in  more  detail.) 
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section  1 

INTRODUCTION  AND  BACKGROUND 


Before  attempting  to  address  the  problem  of  bottom- 
bounce  coherence,  it  is  essential  to  have  a  clear  picture  for 
the  dominant  paths  and  loss  mechanisms  associated  with 
bottom-interacting  signals.  Substantial  progress  has  been 
made  in  this  area  under  the  auspices  of  the  NAVELEX  Bottom 
Interaction  Program  administered  by  NORDA.  In  fact,  the 
understanding  of  loss  mechanisms  has  reached  the  point  that 
an  Advanced  Development,  global,  geo-acoustic  bottom-loss 
model  (BLUG)  is  being  developed  by  SAI  to  support  routine 
Fleet  predictions.  This  paper  attempts  to  summarize  our 
current  understanding  of  the  bottom-loss  problem  and  suggests 
certain  priorities  for  bottom-bounce  coherence. 

The  bottom-loss  problem  has  been  historically 
divided  into  analysis  of  phenomena  where  the  sediments  tend 
to  be  thick  (greater  than  -300  meters)  or  thin  (loss  than 
-100  meters).  in  these  two  categories  bottom  loss  and 
coherence  are  expected  to  have  substantially  different 
characteristics,  in  intermediate-thickness  areas  a  mix  of 
characteristics  is  expected.  Sections  2  and  3  describe  loss 
in  thick  and  thin  sediments,  respectively.  Section  4 
summarizes  on-going  work  at  SAI  in  modeling  the  coherence  of 
these  signals. 

Acknowledgments 

Preparation  of  this  paper  was  supported  by  the 
Naval  Sea  Systems  Command.  The  Bottom  Loss  Upgrade  (BLUG) 
work  was  sponsored  by  the  SEAS  and  TAEAS  Programs,  and  the 
coherence  work  by  the  Bottom  Interaction  Program. 


Section  2 

BOTTOM  LOSS  IN  THICK  SEDIMENT  aREAS 


The  development  o£  a  global  geo^acoustic  model  has 
depended  upon  the  successful  interpretation  of  a  vast  body  of 
bottom“loss  measurements.  Figures  1-4  illustrate  four 
examples  of  bottom-loss  data  collected  by  NADC  in  areas  of 
thick  sediments.  One-octave  bands  were  used  about  the 
indicated  center  frequencies  and  enough  data  were  collected 
th  c  the  curves  represent  median  values  in  2*  bins  in  grazing 
angle.  These  data  are  discussed  in  much  mere  detail  in 
Reference  1. 


These  data  contain  several  consistent  features  and 
a  few  artifacts  of  the  measurement  geometry.  Examples  of 
artifacts  are  the  negative  losses  at  very  low  grazing  angles 
(Figure  1),  and  much  of  the  structure  at  low  frequencies  and 
low  grazing  angles  (Figures  1,  2  and  3).  A  prominent  feature 
in  many  of  the  data  is  the  rapid  increase  in  low-frequency 
loss  around  30-40  degrees  grazing.  This  feature  is  not  an 
artifact,  but  a  mani  station  of  the  fact  that  the  intense, 
shallower  angle  paths  are  not-  reflecting  from  the  water- 
sediment  interface,  but  are,  in  fact,  being  returned  by 
refraction  at  depth  in  the  sediments.  Figure  5  illustrates 
the  ray  geometry  for  this  case.  The  reflected  rays  are 
indicated  by  dashed  rays,  and  inside  the  caustic  range,  R^r 
only  reflected  rays  are  present.  This  leads  to  the  second 
prominent  feature  of  the  data:  the  relatively  angle-  and 
frequency-independent  character  of  the  losses  for  angles 
greater  than  30-40  degrees.  There  is  observed  some  non¬ 
trivial  frequency  dependence  (discussed  in  more  detail  in 
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GRAZING  ANGLE  (DEGREES) 

Rfluf*  1.  Bottom  loaa  veraut  grazing  angle  measured  In  six  one-octave  bands 

at  Indicated  center  frequencies:  zero  high-angle  frequency  dependence. 
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GRAZING  ANGLE  (DEGREES} 

Rgura  2.  Bottom  lots  versus  grazing  angle  measured  in  six  one-octave  bands  at 
indicated  center  frequencies:  por:tive  high-angle  frequency  dependence. 


100  HZ 


(ap)  SS01  vvoiioB 


24-7 


GRAZING  ANGLE  (DEGREES) 


0 

U 


0  0 

c  2 

•>  c 

2  ■ 

O  3 

£  2 
o 


it 

5| 


•I 

£  -o 

E  I 
2  S 

o  I 
CD  .£ 


0 

w 

3 


IZ 


24-8 


GRAZING  ANGLE  (DEGREES) 
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Figiire  5. 


Ray  geometry  (a) ,  angle  (b)  and  intensity  (c) 

versus  range  for  reflected  ( - )  and  refracted 

( - }  paths  where  the  refracted  paths  form  a 

caustic  (*••)  in  the  sediment. 


Reference  1)  reflecting  either  *  pr  ;ence  of  basement, 
and/or  the  effects  of  near-surfa. a  la, erlng  in  the  sedi¬ 
ments  . 


When  one  examines  time  series  of  the  shots  used  to 
collect  .these  data,  this  interpretation  is  reinforced. 
Figure  6  illustrates  the  low  pass  (<2  kHz)  impulse  response 
of  the  returned  wave,  aligned  to  the  first  return  and  plotted 
versus  grazing  angle. ^  The  refracted  arrival  seems  weak 
compared  to  the  initial,  reflected  arrival,  however  hen  this 
time  series  is  low-pass  (40-90  Hz)  filtered  (F  7),  the 
refracted  ar . ival  is  seen  to  dominate.^  The  presence  of  near 
surface  fine-layering  is  seen  in  the  early  structure  in 
Figure  6.  The  constant  spacing  in  time  between  these  micro¬ 
arrivals  suggests  spatially  homogeneous  layers.  In  these 
data  the  corresponding  bottom-loss  measurements  actually  show 
somewhat  less  high-angle  loss  at  high  frequencies  than  at  the 
lower  frequencies  which  essentially  don't  "see*  this  layering 
(Figure  7). 


The  NAOC  data  set  of  several  hundred  bottom  loss 
measurements  were  successfully  interpreted  this  way  in  the 
context  of  a  simplifed  geo-acoustic  model  for  the  sediments 
(Figure  8).  By  examining  various  parts  of  the  bottom-loss 
versus  grazing-angle  curves,  estimates  were  made  of  the 
sound-speed  and  attenuation  profiles,  plus  apparent  surface 
impedances  and  some  aspects  of  the  near-surface  layering. 
Figures  9-12  illustrate  the  comparisons  of  data  (x)  with 
simulated  bottom  loss  (including  estimated  measurement  arti¬ 
facts)  for  four  data  sets.^  These  curves  progress  through 
higher  and  higher  apparent  surface  impedances  (much  higher 
than  the  actual  sound -speed/ density  contrast  at  the  water- 
sediment  interface  would  suggest),  indicating  that  even  at 
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low  fra^ut  icies,  near-surface  layering  nay  be  quite 
important. 


In  summary,  the  behavior  of  loss  in  thick-sediment 
areas  is  dominated  by  the  following  features: 

1.  For  low  frequencies  and  angles  shallower  than 
the  apparent  angle  of  the  caustic  emerging  from 
the  bottom  (30-40  degrees),  the  field  is  con¬ 
trolled  by  the  deep-refracting  path. 

2.  At  higher  frequencies  for  the  same  angles,  the 
more  highly  attenuated  refracting  path  becomes 
dominated  by  the  near-surface  reflected  path, 
at  an  angle  which  depends  on  the  layering 
characteristics. 

3.  At  higher  grazing  angles  the  field  is  con¬ 
trolled  by  the  near-surface  properties, 
including  layering,  and  may  exhibit  a  fairly 
complicated  frequency  dependence  with  only 
slight  angle  dependence. 

4.  As  sediments  get  thinner  (Figure  2)  the  base¬ 
ment  reflected  return  may  dominate  the  low- 
frequency  high-angle  loss.  (For  very  thin 
sediments  this  return  may  dominate  even  up  to  1 
kHz  as  discussed  in  the  following  subsection.) 
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The  same  bottom-loss  simulation  model  has  been  run 
against  bottom-loss  data  collected  in  several  areas  of  the 
Pacific  where  the  sediments  tend  to  be  approximately  40 
meters  thick.  Figures  13-15  illustrate  some  examples  of 
these  comparisons  which  tend  to  show  less  agreement  for 
several  reasons  discussed  below. ^  The  data  generally  show 
less  loss  at  lower  frequencies  and  very  little  obvious  angle 
dependence/  except  at  very  shallow  grazing  angles  where  the 
path  can  refract  above  the  basement. 
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In  these  thin-sediroent  areas  the  basement  is  quite 
rough  and  large-angle  scattering  is  expected.  Figure  16 
illustrates  the  aligned  low-pass  (<1  kHz)  impulse  response 
from  such  an  area.*^  Note  the  time-axis  scale  change  over 
Figures  6  and  7.  The  signal  is  reverberant  from  the  bottom/ 
lasting  well  over  a  second.  At  low  frequencies  (Figure  17) 
the  signal  lasts  just  as  long  with  no  dominant  path.*^  The 
scattering  must  be  quite  diffuse/  since  the  "specular” 
(first)  return  contains  very  little  energy. 

The  fact  that  the  time  spreads  at  low  and  high  fre¬ 
quencies  are  similar  suggests  that  these  returns  cannot  be 
coming  from  great  depths  in  the  basement  itself  (or  no  high- 
frequency  returns  would  be  seen).  Hence  the  signal  must  be 
coming  from  the  sediment/basement  interface  (energy  analysis 
shows  it  cannot  be  a  pure  sediment-reflected  return).  When 
the  observed  time  spread  is  translated  into  the  corresponding 
area  on  the  bottom  from  which  the  signal  must  be  coming 
(Figure  18)/  the  constant  time-spread  ellipses  covering  1  to 
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1*5  seconds  ace  quite  large.  Hence  the  bo ttom-'ce flecked 
signal  is  spread  in  both  azimuthal  and  depression  angles  aa 
well  as  time.  I£  the  source  and/or  receiver  of  Figure  18 
were  moving /  these  angular  spreads  would  result  in  a  fre** 
quency  spread  for  a  cw  signal  due  to  differential  Doppler  on 
the  angularly  spread  components.  The  following  section 
discusses  these  coherence  issues  in  more  detail. 
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Section  4 

COHERENCE  MODELING 

Figure  19  summarizes  the  three  bottom-interaction 
coherence  problems  SAI  is  currently  pursuing  in  the  Bottom 
Interaction  Program.  In  thick-sediment  areas  two  issues  are 
being  addressed: 

1.  The  loss  and  toraporal  spread  of  signals 
reflected  by  the  near-surface  layering.  Greene 
has  just  begun  work  on  this  problem.  In  his 
model  the  layering  is  treated  as  a  discrete  ^ 
jump  (Poisson)  process  with  very  large 
(effectively  infinite)  horizontal  correlation 
lengths.  ^  Besieris  and  Kohler^  are  treating  the 
same  problem  with  continuous  layering  (which 
may  be  slightly  less  physical),  but  with  a 
finite  horizontal  correlation  length. 

2.  The  spatial  coherence  (angular  spreads)  of  the 
refracting  path  associated  with  random 
inhomogeneities  in  the  bulk  of  the  sediment. 
(Dozier  discusses  this  problem  in  more  detail 
in  a  subsequent  paper.) 

In  thin-sediraent  areas  the  loss  and  coherence 
problems  are  not  so  easily  separated.  Our  work  in  this  area 
addresses  all  spreads  as  well  as  total  loss.  The  approach 
uses  a  transport-equation  model  to  construct  the  average 
intensity  as  a  function  of  angles  (0,4)  and  time  (t)  for  a 
given  range  (R)  and  frequency  (f): 


THICK  SEDIMENTS 


THIN  SEDIMENTS  -  ROUGH  BASEMENT 


Transforms  of  this  spread  function  yield  the  spatial  and 
frequency  mutual  coherence  functions. 

In  terms  of  observables  and  sonar  equation 
parameters,  for  a  cw  signal  (static  source  and  receiver) 

“  /  l(6/^,t)  dt 
t 

defines  the  angul.ir  spectrum  of  the  signal  and  a  corre¬ 
sponding  beam  output  for  steered  direction  6s»^s  would 
be 

^B^®s'*s)  “  ^  ^  BP(0,*;0  ,4  )  d0d^ 

0  ^  as 

where  BP  is  the  beam  pattern  for  the  steered  direction.-  (For 
a  moving  source/receiver  the  corr  ^spending  frequency  spread 
would  be  given  by  the  differential  Doppler  on  these  spread 
angles  computed  at  both  ends,  if  necessary,) 

To  compare  with  impulse  response  data  similar  to 
Figures  16  and  17  we  integrate  over  angle 

I^(t;R,f)  »  /  /  I{6,*,t;R,f)  d0d^ 

e  ^ 

Because  scattering  is  so  dominant,  any  model  for  the  total 
loss  must  be  consistent  vith  these  observed  spreads. 
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The  computation  of  this  basic  intensity  spread 
(l(0,^,t))  is  being  done  by  a  brute-force  approach  in  which 
the  bottom. is  divided  into  small  areas  related  to  the  source 
and  receiver  by  two  grazing  angles  (0,6')  and  a  bistatic 
scattering  angle  (A^).  (The  bottom  roughness  spectrum  is 
taken  to  be  isotropic,  hence  only  A*  is  important.)  These 
paths  are  related  through  refraction  in  the  sediment  to  the 
corresponding  basement  angles  (0jj#6^)#  and  appropriate  (fre¬ 
quency-dependent)  attenuation  losses  in  the  sediment  are 
accumulated. 

The  power  scattered  from  into  6^-  A#  is  computed 

using  a  composite-surface  scattering  kernel  developed  by 
Brown  (Figure  20).®  The  composite-surface  theories  recognize 
that  the  large-angle  scattering  from  a  surface  is  controlled 
by  the  portion  of  the  surface  spectrum  where  the  acoustical 
(ko)  and  surface  wavenumbers  are  comparable  (basically 
diffractive  Bragg-like  scatter).  These  "diffraction* 
gratings  are  essentially  carried  on  the  large  (long- 
wavelength)  surface  which  introduces  tilts  to  the  gratings 
and,  at  high-frequencies,  can  shadow  parts  of  the  surface. 
The  analysis  is  done  in  terms  of  these  two  surfaces  which  are 
separated  spectrally  at  k*,  which  depends  on  kg.  The 
results  are  fairly  insensitive  to  the  precise  value  of  k*. 

The  corresponding  scattering  kernel  then  transi¬ 
tions  smoothly  between  the  physical  optics  (high-frequency) 
and  diffractive  (low-frequency)  limits.  We  are  in  the 
process  of  applying  this  model  to  analyze  several  data  sets 
both  for  loss  and  measured  spreads. 
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Figure  20.  Composite-surface  scattering  from  basement. 
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ABSTRACT 


Deep  Tow  data  from  an  abyssal-hill 
site  in  the  Northeastern  Pacific  (north 
of  Hawaii)  have  been  analyzed  to  estimate 
the  statistics  of  basement  rouahness  and 
sediment  thickness.  Specifically/  base¬ 
ment  rouqhness  spectra  have  been  devel¬ 
oped,  tested  for  directionality  and 
stationarity,  and  compared  with  Bell's 
results  [1975]  for  a  similar  area  but 
larqer  wavelenqths.  Distribution  func¬ 
tions  of  basement  heiqht  (detrended)  and 
slope  and  sediment  thickness  are  also 
presented. 

Ocean  bottom  rouqhness  is  a  quantity  of  interest  to  both 
the  geophysical  and  acoustical  communities.  Rouqhness  varies 
from  large  scale  features  which  are  essentially  bathymetry  to 
features  whose  scale  lengths  are  so  short  that  only  statisti¬ 
cal  descriptions  are  appropriate.  Geologists  have  found 
roughness  a  useful  description  of  the  ocean  bottom.  In 
addition,  oceanographers  have  shown  an  interest  in  rouahness 


as  a  generator  of  open-ocean  Internal  waves  [Bell,  1975],  and 
acousticians  have  seen  roughness  as  a  mechanism  for 
scattering  of  bottom-incident  sound  waves  [Berkson  and 
Matthews,  1984]. 

In  the  abyssal  hills  of  the  North  Pacific  at  acoustical 
frequencies  from  tens  to  thousands  of  hertz  the  bottom  prop¬ 
erties  of  most  importance  are  the  sediment  thickness  to  the 
first  major  internal  reflector  and  the  roughness  of  that 
interface.  The  relatively  small  impedance  contrast  at  the 
water-sediment  interface  produces  low  reflectivity  compared 
to  that  caused  by  the  large  contrast  at  this  strong  reflector 
in  this  frequency  regime.  The  high  impedance  contrast  at 
this  reflecting  horizon  {acoustic  basement),  together  with 
its  roughness,  tends  to  scatter  back  into  the  water  column  a 
large  amount  of  energy  over  relatively  wide  angles  [Spofford, 
1983;  Diachok  et  al.,  1984]. 

The  acoustical  properties  of  interest  for  this  problem 
are  the  various  wave  speeds  and  attenuations,  the  impedance 
contrasts,  the  sediment  thickness,  and  the  scatter¬ 
controlling  aspects  of  basement  roughness.  This  paper  is 
concerned  with  estimates  of  the  last  two  properties  for  an 
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abyssal-hill  environment  representative  of  larae  reqions  of 
the  North  Pacific.  The  data  used  to  address  these  issues 
were  collected  durinq  a  Deep  Tow  survey  fSpiess  and  Lonsdale, 
1982]  of  a  manqanese  nodule  site  in  the  North  Pacific  known 
as  MANOP  Site  R.  The  oriqinal  qoals  of  the  expedition  were 
to  map  the  local  variability  of  the  extensive  nodule  coveraqe 
at  this  latitude  usinq  near-bottom  side-lookinq  sonar  and 
photography,  and  to  determine  sites  which  would  be  appro¬ 
priate  for  subsequent  deployment  of  bottom  landinq  vehicles 
to  carry  out  experiments  in  the  sediment.  Site  R  was 
selected  as  typical  of  the  extensive  nodule-covered  red  clay 
region  of  the  North  Pacific.  The  total  sediment  column,  as 
determined  by  seismic  refraction  stations  within  a  few 
hundred  miles,  ranqes  from  300  to  500  m  [S.  M.  Smith,  SIO, 
personel  communication).  The  first  maior  internal  reflector, 
constituting  the  effective  "basement"  for  the  acoustic 
aspects  considered  here,  lies  at  a  depth  of  30  to  50  m  below 
the  sea  floor.  The  dense  surficial  manqanese  nodule  coveraqe 
makes  the  sea  floor  itself  the  maior  reflector  at  frequencies 
above  a  few  thousand  hertz,  and  at  4  kHz  there  are  two 
additional  weak  reflectors  discernable  above  the  effective 
acoustic  basement  (Spiess  and  Weydert,  1984). 
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Data  from  over  100  km  of  tracks  crisscrossing  a  200  sq 
km  area  were  analyzed  with  the  primary  objective  of 
describing  the  topographic  roughness  of  the  major  internal 
reflector  and  a  secondary  objective  of  characterizing  the 
thickness  of  the  sediment  overlying  it. 

Most  data  on  seafloor  roughness  are  collected  with 
various  instruments  at  (or  very  neat)  the  ocean  surface  and 
have  limited  spatial  resolution.  Berkson  and  Matthews  [1984] 
summarize  data  that  typically  resolve  features  with  scale 
lengths  of  200  to  2000  m.  Only  scale  lengths  greater  than 
half  the  acoustical  wavelength  can  contribute  to  the  scat¬ 
tered  field;  hence  most  of  the  previous  roughness  data  are 
not  applicable  to  the  basic  scattering  problem.  However, 
when  the  roughness  power  spectrum  appears  to  follow  a  simple 
power  law,  extrapolation  to  muc  shorter  wavelengths  might 
provide  the  necessary  data. 

The  roughness  issuer  that  this  analysis  specifically 
addresses  are: 


(1)  Pow^er  spectra  at  wavenumbers  of  acoustical  interest 
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(2)  Directionality  of  basement  spectra 

(3)  Spatial  stationarity 

(4)  Consistency  with  longer  wavelength  data  (including 
the  extrapolation  issue) 

(5)  Height  distributions 

(6)  Slope  distributions. 

In  addition#  distributions  of  sediment  thickness  have  been 
computed. 

The  following  sections  describe#  respectively#  the 
measurement  apparatus,  the  data  processing  and  analysis,  the 
results#  and  the  conclusions  regarding  these  issues. 


MEASUREMENT  APPARATUS 

The  Deep  Tow  Instrumentation  System,  originally  deployed 
in  1962#  has  undergone  a  series  of  augmeiuations#  allowing  it 
to  examine  with  increasing  resolution  and  diversity  the 
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properties  of  the  sea  floor.  Typically  Deep  Tow  (fish)  is 
towed  10-100  m  above  the  sea  floor.  A  12-kHz  transponder  net 
is  deployed  on  the  bottom  allowing  the  tow-ship  and  the  fish 
to  determine  their  positions  relative  to  the  bottom  with 
accuracies  of  20  m  and  2  m  respectively. 

The  recordings  of  interest  to  this  analysis  are  made 
using  the  4-kHz  seismic  profiler  [details  may  be  found  in 
Lonsdale,  Tyce  and  Spiess,  1974],  Although  this  system  is 
effectively  omnidirectional,  it  is  so  close  to  the  bottom 
that  it  has  an  effective  resolution  of  a  few  meters  for  the 
roughness  seen  in  this  area. 

The  soundings  are  virtually  continuous  and  produce  a 
graphic-recorder  depth  chart  for  the  fish  and  for  all  reflec¬ 
tors  (in  this  case  the  water-sediment  and  sediment-basement 
interfaces)  versus  time.  From  the  geometry  reconstructed 
using  the  cransponders,  time  is  converted  to  position. 

DATA  PROCESSING  AND  ANALYSIS 

Figure  1  illustrates  the  entire  Deep  Tow  track  at  the 
Manop  "R"  site.  The  origin  (0,0)  of  the  local  coordinate 


system  was  30*19. 9'N,  157*50. 9'W.  The  heavy  dashed  and  solid 
portions  of  the  track  were  digitized.  The  portions  not 
digitized  had  poor  trace  quality  or  too  many  interruptions  to 
be  useful  for  spectral  estimates.  The  portions  digitized  but 
not  analyzed  had  poor-quality  information  on  the  fish  depth. 
The  numbering  scheme  corresponds  to  specific  analysis 
comments  which  follow. 

For  data-processing  purposes,  the  fish  track  segments 
were  smooth  enough  that  we  assumed  h{r)  =  h(s)  where  r  is  a 
straight-line  range,  s  is  path  length  along  the  actual  track, 
and  h  is  depth.  Probability  density  functions  of  slope  and 
sediment  thickness  were  computed  using  unfiltered  data. 
Basement  depth  distributions  were  computed  after  detrending, 
effectively  simulating  a  high-pass  filter  (Figure  2).  Note 
that  the  sea  floor  depth  profile  appears,  in  a  sense,  to  be  a 
low- pass  filtered  version  of  the  topography  of  the  first 
major  reflector. 

In  spectral  analysis,  the  high-pass  wavenumber  corre¬ 
sponded  to  a  wavelength  approximately  one-tenth  the  length 
of  the  track  being  processed.  Low-pass  wavenumbers  ranged 
from  .008  m"^  for  Track  8  to  .025  m~^  for  Track  6.  One¬ 
dimensional  power  spectral  densities  were  computed  using  an 
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FFT  from  windowed,  detrended  data.  The  window  consisted  of  a 
raised  cosine  applied  to  the  first  and  last  one-eiqhth  of  the 
points  on  each  track.  The  soectrum  of  the  window  has  a  side- 
lobe  rolloff  (P(K)  ~  much  steeper  than  the  spectra 
of  the  data,  and  a  narrower  main  lobe  than  a  Hanninq  window. 
Transform  size  depended  on  track  lenqth.  Ranqe  steps  were 
approximately  3.5  meters,  orovidinq  spectral  estimates  to  a 
maximum  wavenumber  of  about  1  m"^. 

Where  tracks  were  lonq  enouqh,  averaqes  of  adiacent 
segments  of  a  track  were  used  to  produce  a  smoothed  spectrum. 
In  these  cases,  a  runs  test  was  used  to  determine  a  minimum 
subtrack  lenqth  so  that  adjacent  subtracks  were  uncorre¬ 
lated. 


RESULTS 

Track  8  (see  Figures  1  and  2)  was  chosen  to  illustrate 
results  from  individual  tracks.  The  oower  soectral  density 
for  this  track  (Figure  3)  can  be  aporoximated  bv  P(K) 
K“*^,  where  P  is  the  one-dimensional  power  soectral  density, 
K  IS  wavenumber  and  b  characterizes  the  slope  (on  loq  axes) 
of  the  h  gher  frequency  portion  of  the  spectrum.  For  Track 


8,  b=3.8.  Probability  density  functions  of  detrended  base¬ 
ment  depth,  one-dimensional  basement  slope,  and  sediment 
thickness  are  shown  in  Fiqure  4.  The  one-dimensional  rms 
basement  slope  for  this  track  is  4.8  decrees.  A  qaussian 
slope  distribution  would  imply  an  equivalent  two-dimensional 
rms  slope  of  6.8  deqrees. 

Table  1  summari25es,  for  all  tracks,  the  one-dimensional 
slopes,  rms  rouqhnesses,  and  values  of  b.  Values  of  b  were 
nearly  identical  for  all  the  tracks  (Fiqure  5)  whereas  rms 
rouqhness  and  rms  slope  varied  by  almost  a  factor  of  2  (Table 
1).  Clearly  over  this  area  these  statistical  measures  are 
not  constant,  indicatinq  some  deqree  of  non-stationaritv  of 
the  random  process  h(x)  which  miqht  be  relevant  to  acoustical 
scatterinq  effects.  Tracks  2a,  6a  and  8a  (see  Fiqure  1)  were 
used  to  test  for  directionality  of  basement  orooerties. 
Fiqure  6  illustrates  the  correspondinq  power  spectra.  Within 
the  limits  imposed  by  the  relatively  small  samples  it  is 
difficult  to  discern  a  siqnificant  directionality. 

Assuminq  approximate  statistical  homoqeneity  over  the 
area,  an  enserable-averaqe  spectrum  was  computed  (Fiqure  7). 
In  this  region  it  is  clear  from  the  individual  and  ensemble 
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spectra  that  in  the  wavenumber  domain  from  iqO  jh'I 
the  spectral  slope  steepens  substantially.  This  is  not  an 
artifact  of  the  data  processing  or  the  omni-directional 
nature  of  the  transducer.  (The  slight  decrease  in  slope  near” 
10^  m“^  might  be  a  sampling  artifact,  or  it  might  represent 
the  beginning  of  a  new  transition  to  fine-scale  roughness.) 

In  Berkson  and  Matthews  [1984]  other  abyssal-hill  areas 
of  the  Pacific  show  power  spectra  with  slopes,  b,  between  2.2 
and  4.9  for  wavenumbers  of  .003  to  .03  m*^.  Bell  [1975] 
reported  a  spectrum  representative  of  over  8000  km  of  pro¬ 
files  from  published  surface-source  data  and  data  taken  by 
Mudie  using  Deep  Tow  in  a  different  abyssal-hill  area. 
These  data  have  been  plotted  in  Figure  8  along  with  our 
ensemble  spectrum.  Both  the  slopes  and  spectral  levels  in 
the  common  regions  are  in  excellent  agreement.  Figure  8 
clearly  illustrates  the  danger  in  extrapolating  the  shorter 
wavenumber  spectra  to  the  higher  wavenumbers  of  acoustical 
interest  using  only  the  spectral  slope  of  the  short- 
wavenumber  data.  However,  the  spectra  derived  from  Bell's 
data  (from  a  wide-area  average)  and  from  the  data  of  the 
present  study  (from  a  relatively  small  area)  match  closely  in 
the  wavenumber  region  of  overlap  (.006-, 03  ra"^).  This  match 
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suggests  that,  for  wavenumbers  greater  than  .006  the 
basement  roughness  spectrum  may  be  relatively  universal.  The 
agreement  also  suggests  that  below  this  wavenumber  the  true 
basement  roughness  and  bathymetry  have  comparable  roughness. 

Finally,  Figure  9  displays  ensemble  distribution  func¬ 
tions  of  basement  height,  basement  slope,  and  sediment  thick¬ 
ness.  Overplotted  on  the  histograms  are  Gaussian  distribu¬ 
tions  with  the  same  means  and  variances  as  the  corresponding 
data.  The  tendency  for  all  three  Gaussian  curves  to  be 
broader  than  the  apparent  distribution  functions  means  that 
the  tails  of  the  histograms  are  driving  the  variance  of  the 
Gaussians.  For  slope  and  sediment  thickness  these  differ¬ 
ences  are  not  too  important.  The  clearly  non-Gaussian  height 
distribution  can  pose  significant  problems  for  acoustical 
scattering  theories,  nearly  all  of  which  assume  a  Gaussian 
distribution  for  h(x). 

SUMMARY  AND  CONCLUSIONS 


Deep  Tow  data  from  an  abyssal  hill  area  of  the  Pacific 
north  of  Hawaii  have  been  analyzed  to  provide  statistics  of 
basement  roughness  and  sediment  thickness.  Detrended 
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basement  roughness  shows  rms  height  variations  of  between  2 
and  6.5  m  depending  on  the  specific  track  in  the  region 
surveyed.  Basement  rms  slopes  varied  between  3.7  and  8.6 
degrees  with  tougher  basements  tendinq  to  show  larger  slopes. 
The  height  distribution  is  distinctly  non-Gaussian / 
suggesting  significant  consequences  for  bottom-bounce 
scattering  theories.  The  slope  distribution  is  more  nearly 
Gaussian. 

Basement  roughness  spectra  for  the  wavenumber  domain 
10“2  to  10°  m~^  show  a  range  of  logarithmic  slopes  between 
-3.5  and  -4.0.  Little  directionality  was  found  in  the  power 
spectra  for  a  set  of  three  crossinq  tracks.  In  the  wave- 
number  region  of  data  overlap,  there  was  significant  agree¬ 
ment  in  both  level  and  spectral  slope  (-2.5)  between  the 
ensemble-averaged  power  spectrum  from  the  Deep  Tow  data  and 
Bell's  (19751  large-scale  Pacific  power  spectrum.  The 
present  study  extends  the  basement  power-spectrum  measure¬ 
ments  upwards  by  two  wavenumber  decades.  These  results 
should  be  applicable  to  acoustical  scattering  problems  for 
frequencies  up  to  at  least  a  few  hundred  hertz. 
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Table  1.  Track  and  ensemble  summary  statistics 
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Sample  track-data  for  Track  8  of  Figure  \ 
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Power  spectral  densities  of  detrended  basement 
depth  for  all  8  tracks. 


Fig.  6.  Power  spectral  densities  of  crossing  tracks  to 
test  potential  directionality. 
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Fig.  8.  Comparison  of  power  spectral  densities  for 
low  wave-numbers  (Bell,  1975)  and  composite 
spectrum. 


Fig.  9.  Ensemble  data  histograms  for  all  tracks  with 
corresponding  Gaussian  fits;  (a)  detrended 
basement  depth;  (b)  basement  slope;  and 
(c)  sediment  thickness. 
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